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Making use of recently released data on dihadron correlations by the STAR collaboration, I analyze 
the long-range ( "ridge-like" ) part of these data and show that the dependence on both transverse 
momentum as well as orientation with respect to the event plane are consistent with correlations 
expected from only collective flow. In combination with previously analyzed centrality-dependent 
data, they provide strong evidence that only collective flow effects are present at large relative 
pseudorapidy. In contrast, by analyzing a "background subtracted" signal, the authors presenting 
the new data concluded that the ridge-like part of the measured correlation could not in fact be 
entirely generated from collective flow of the medium. I explain the discrepancy and illustrate some 
pitfalls of using the ZYAM prescription to remove flow background. 



I. INTRODUCTION 

Two-particle correlations in relativistic heavy ion colli- 
sions — the probability of seeing a pair of particles with 
relative azimuth Acf> = (0i — <f>2) and relative pseudo- 
rapidity A77 = {rji — 772) in a given collision event - 
display unique features not seen in other types of colli- 
sions such as p-p or d-Au. In particular are significant 
long-range correlations extending to large A77, which of- 
ten show interesting ridge and shoulder structures after 
a model-dependent subtraction of elliptic flow [J-Q- 

Recently it was proposed that this long-range part of 
the correlation could be entirely explained by collective 
flow effects [3, Central to this idea was the fact 
that, due to event- by- event fluctuations, there should ex- 
ist not only elliptic flow, but also "triangular flow" , which 
should add a non-negligible contribution to these data. 
Although triangular flow (w 3 ) has not yet been directly 
measured at RHIC, there were hints from a transport 
model, and it was later shown Q from viscous hydrody- 
namic calculations, that the centrality dependence as well 
as the size (depending on viscosity) of the third Fourier 
component of the correlation V3A = (cos(3A0)} does in- 
deed quantitatively match that expected as arising sim- 
ply from triangular flow. Thus, the measured dihadron 
correlation at large Aij, consists almost entirely of the 
lowest few Fourier components, each of which can be 
quantitatively understood as coming from collective flow 
(plus global momentum conservation) — at least for the 
data analyzed, which had transverse momentum triggers 
of p t = 2.5 GeV and lower. 

Even more recently, the STAR collaboration has re- 
leased data from 200 A*GeV Au-Au collisions at RHIC, 
where a trigger particle (with a p t of 3-4 GeV or 4-6 
GeV) was restricted to be at a fixed angle with respect 
to the measured event plane, and its correlation with 
associated hadrons in various pt bins was measured Q. 
The data were then separated into "ridge-like" and "jet- 
like" components, where the ridge-like correlation was 
defined by a projection on |A?7| > 0.7, while the jet-like 
correlation was defined by taking the total short-range 
correlation at |Ary| < 0.7 and subtracting the ridge- like 
correlation. 



These new data provide an opportunity to test the idea 
that the long-range ridge-like correlation may be gener- 
ated exclusively by collective flow effects. 

In Ref. 0] , the ridge- like correlation was analyzed using 
the fairly common Zero- Yield- At-Minimum (ZYAM) pre- 
scription [8[ to "subtract" the elliptic flow signal which 
dominates the unsubtracted data, in order to study non- 
flow behavior of the system. These flow-subtracted data 
show a dependence on the angle between the trigger par- 
ticle and the event plane, and this was taken as evi- 
dence that collective flow alone could not explain this 
signal since triangular flow should be uncorrelated with 
the event plane. 

In this article I show that this conclusion is not cor- 
rect. As was seen in previously analyzed data, the long- 
range "ridge-like" correlation is entirely consistent with 
what is expected purely from collective flow, and rather 
than contradicting this idea, the new data actually pro- 
vide strong evidence. I then illustrate how the misleading 
background-subtracted signal seen by STAR can be gen- 
erated by the use of the ZYAM prescription in a situation 
where the assumptions of such a subtraction procedure 
do not hold. 



II. FOURIER DECOMPOSITION OF 
LONG-RANGE CORRELATIONS 

It is useful to do a Fourier decomposition of unsub- 
tracted dihadron correlation data at large relative pseu- 
dorapidity. These data contain only long-range correla- 
tions, which stem from early times in the collision evolu- 
tion [H , and it has been argued [|[ that they may contain 
significant contributions only from the collective behav- 
ior of the system (along with a trivial correlation from 
global momentum conservation). The recent data from 
STAR present an opportunity to further test this idea by 
adding information about orientation with respect to the 
event plane as well as p% dependence, which complements 
the previous analyses 0, |g| of centrality dependent data 
from PHOBOS 0, [Io| and STAR [n|. 

A Fourier analysis was done of the ridge-like compo- 
nent of the recent data (i.e., correlations for \Ar/\ > 0.7) 
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FIG. 1. (Color online) First 4 harmonics V„a = (cos(nA^))) 
of dihadron correlations for particle pairs with relative angle 
Acj> — 4> t — 4>a, and relative pseudorapidity |A?7| > 0.7 Q. On 
the left VnA is plotted as a function of angle <f> s between a 3-4 
GeV trigger particle and the event plane, and on the right as 
a function of transverse momentum of the associated particle, 
p[ a ^ . The even harmonics also include bands for the estimated 
flow background [jj. Points are placed at the midpoint of each 
bin. The grey curve of V3A versus pt shows the p t dependence 
of triangular flow from viscous hydrodynamics (see text for 
details). Data for 4-6 GeV trigger particles show the same 
trends, but with larger statistical fluctuations. 



Q. At this 20-60% centrality range, the correlation 
is dominated by the second Fourier component V2A = 
(cos(2A0)) (unless the angle of the trigger particle with 
respect to the event plane is close to 7r/4). This is fol- 
lowed by the remainder of the first four harmonics (V4A, 
V3A, and Via) while all higher moments are significantly 
smaller — the data in every bin can be very well approx- 
imated by the first four Fourier coefficients, 



dN 



N 



dA(j) 2tt 



[1 + 2Via cos(A^) + 2V" 2 a cos(2A0) 
+ 2V3A cos(3A</>) + 2 V 4A cos(4A0)] . (1) 



In general, collective flow introduces no direct correla- 
tion between pairs of particles. In the absence of non-flow 
correlations, the two-particle correlation is determined 
solely by the one-particle distribution, and so in a given 
collision event the flow contribution to each harmonic 
factorizes: 

V nA = (cos(nA</>)> = Re(e m( ^-* t) ) 
=Re (y»^a-^)^ e -™W t -^)^ 

=(cos(n4> a - nipn)) (cos {nipt - nip n )) (2) 

Here 4> a and (fit are the angles of the associated and 
trigger particle, respectively. ip n is an angle that can 
be determined separately for each n in a given event in 
the same way as the event plane angle for elliptic flow, 
1P2 = ipEP- Due to event-by-event fluctuations, these 
angles are not in general equal to each other or to the 
reaction plane that is defined by the beam and impact 
parameter. Note also that when the correlation is av- 
eraged over many events, the average of the product is 
sensitive to the magnitude of these fluctuations. 

Specifically how these two factors depend on indepen- 
dently measured quantities (e.g., V2, etc. depends on how 
the pairs are selected and averaged. This is discussed for 
each individual Fourier harmonic in the following sec- 
tions: 



The even harmonics: elliptic and quadrangular 
flow 



The dominant harmonic, V2A, is understood as coming 
largely from elliptic flow. The flow contribution is 12]: 



V 2A =(cos(2A(/))) = v. 



(3) 



Here v. 



(cos(20 a — 24>ep)} is the elliptic flow of par- 



ticles in the same bin as the associated particle, while 
y {t,R) ^ cos (2^ s ) since the trigger particle is fixed to be 
at a particular angle 4> s = \<j>t — iPep\ with respect to the 
event plane, but there are corrections due to the finite 
angular bin width and event plane resolution [l3| . These 
corrections depend on vi and 1)4 of particles in the same 
p t bin as the trigger particle (see Eq. (4) of Ref. 0] for 

(t Ft) 

the precise functional form of v 2 used by STAR) . 
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Figure [Tj (row 2) shows the extracted V2A along with 

the values of v^'v^'^ reported by STAR (with bands 
representing their estimated systematic uncertainty). By 
comparing the extracted values and estimated back- 
ground, one can see that the data show the behavior 
expected from collective flow both as a function of pt 
and trigger angle. The only difference is that the mag- 
nitude of the estimated flow contribution is consistently 
smaller than the extracted Fourier component — the im- 
plied non-flow signal has roughly the same pt dependence 
as flow, and contributes a correlation that varies mono- 
tonically as the orientation goes from in-plane to out-of- 
plane; i.e., it has the same behavior as elliptic flow. 

So either the background estimation is correct and the 
non-flow contribution to (cos2(A(/>)) has the same prop- 
erties as flow, or the flow background is underestimated. 
In fact, the latter is quite likely when the method of ob- 
taining the estimation is analyzed. 

The lower bound of the V2 systematic uncertainty band 
for both the associated and trigger particle are given by 
measurements of the four particle cumulant elliptic flow 
«2{4}. This measurement significantly reduces the effects 
of non-flow [l4| , but also has a negative contribution from 
elliptic flow fluctuations (l5j . Because of the imperfect 
event plane resolution, however, the correlation V2A here 
has a positive contribution from elliptic flow fluctuations, 
and so the lower edge of the background uncertainty band 
is certainly smaller than the actual flow contribution. 

The upper bound of the uncertainty band is obtained 
by using an event plane V2{EP} for the trigger particle 
and an away-side two-particle cumulant i>2{2, AS} for the 
associated particle, though the dominant contribution is 
from the latter (recall that vi of the trigger particle only 
comes in as a small correction due to finite bin width and 
event plane resolution). 

The away-side two-particle cumulant is obtained in 
essentially the same way as one would obtain the sec- 
ond Fourier component of a two-particle correlation, ex- 
cept the correlation is only integrated over the away side 
(|A</>| > instead of the entire azimuth. This is pre- 
sumably done to remove the non-flow contribution of the 
jet peak at small A(f> (and small A77). However, as a con- 
sequence, one is no longer actually picking out the second 
harmonic. As mentioned in Ref. [7[, the negative dipole 
Via tends to increase i>2{2, AS}. However, as can be seen 
by integrating Eq.[Tjover \A<j>\ > a third harmonic V3A 
will decrease the value by an even larger amount (as long 
as it is at least as large as §ViA) which is expected and 
also coincides with observation). 

It is quite likely, then, that even the upper uncertainty 
bound of the estimated background is smaller than the 
actual elliptic flow contribution. 

According to Ref 0, V4A has a similar contribution 
from collective flow: 

F 4A = (cos(4A^)~4 o) 4* ,il) , (4) 

where — (cos(40 a — 4:ijj EP )) is the associated quad- 
rangular flow measured with respect to the event plane 



and Vj^ R ^ behaves roughly as cos(40 s ) (see Eq. (5) of 
Ref. Q). Actually, this assumes that in each event 
the 4th harmonic of the single particle distribution is 
centered about the same angle as the 2nd harmonic 
(ipi = tpEp)- In reality, the particle distributions de- 
fine an angle tp4, that, like the event plane, fluctuates 
around the reaction plane event- by- event. The actual 
contribution can be written as Eq. Q plus a somewhat 
complicated additive correction that depends on event- 
by-event fluctuations. So far, U4 has only been measured 
with respect to the event plane, and not with respect to 
its own angle tp4, and a more complete analysis is left to 
future work, but here it is important to note only that 
this correction is expected to be relatively small, but al- 
ways positive. 

Figure Q] (row 4) shows V4A as well as the expected 

background v^v^' R ' as used in Ref. Q. One can see by 
comparing the extracted curves to the bands representing 
the expected background that it behaves as expected as 
a function of p t and angle, and actually shows evidence 
of ip^ fluctuations, since the extracted curves are always 
slightly above the estimated background. Any non-flow 
contribution would need to have the same angular and 
Pt dependence as flow to see these observed trends, or 
it would need to be very small and mimic the effects of 
■04 fluctuations. Note that this means any non-flow sig- 
nal must have an angular dependence for the 4th Fourier 
harmonic that is quite different from that of its 2nd har- 
monic. 

Thus, the observed even harmonics are consistent with 
those expected from flow alone, according to every pos- 
sible test criterion. 



B. The odd harmonics: triangular flow, momentum 
conservation, and directed flow 

The third component, V3A, will have a contribution 
from collective flow of the form 0, @| 

V 3A = (cos(3A0)) ~ vg>v$\ (5) 

where U3 — (cos(30 a — Sips)) is the triangular flow of 

the associated particle, and Ug is similarly the triangular 
flow for the trigger particle. These are flow parameters 
that could be measured with respect to an event- by- event 
plane of triangularity ^3 that is expected to be largely 
uncorrelated with the event plane 0, [la, LLZl , and so this 
contribution should not depend on the event plane angle. 

Indeed, the data show no systematic dependence on 
event plane orientation (see Fig. [TJ row 3). In addition, 
although no direct measurement of 1*3 has yet been made, 
both the size and pt dependence are consistent with hy- 
drodynamic predictions (the solid grey curve in Fig. Q] 
shows the ^-differential v% for a viscous hydrodynamic 

calculation of a 30-35% central collision with n/s = 4~ 

— 1 11 Jin 

and Glauber initial conditions from Ref. [6], multiplied 
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by 0.09, representing the estimated V3 of the trigger par- 
ticle at 3-4 GeV). Note that the centrality dependence 
of this correlation has also been shown to follow that 
calculated in hydrodynamic models @. 

Thus, if there is a non-flow signal, it must have a 
(cos3(A0)} component that is independent of angle, and 
has a p t dependence that is not too different from trian- 
gular flow. 

The first harmonic, Via, has a contribution from global 
momentum conservation shared between a finite multi- 
plicity 18]. Although this is not, strictly speaking, a flow 
correlation, it will be present in any finite-multiplicity 
system, and so is not an interesting non-flow signal. It 
does not have any dependence on the angle with respect 
to the event plane but should be negative with a linear 
dependence on p t of the associated particle [HI . 

In addition, there is a possible contribution from col- 
lective flow. The directed flow that has been measured at 
RHIC with respect to the impact parameter is very small 
everywhere within the acceptance of the STAR TPC (l9j . 
Thus, one might expect no contribution from collective 
flow. However, there can also be a v\ that is generated 
event- by- event due to fluctuations in the initial geome- 
try [20( ] . This contribution is expected to be approxi- 
mately independent of rapidity, and therefore cancels in 
the existing directed flow measurements, which are odd 
in rapidity by construction [2l|,[22[. However, it does not 
cancel in this dihadron correlation, and the effect is given 
by [HHI 

Via = (cos(A0)) ~v[ a) v[*\ (6) 

where V\ = (cos(</> — ipi)) is the directed flow that can 
be defined with respect to an event-by-event angle ipi- 
Generically it should be negative at low p t and positive at 
high p t [20|, with a zero crossing at ^800 MeV such that 
the net momentum in the direction of ipi is approximately 
zero. As mentioned, it is not directly correlated with 
the event plane, so, like the correlation from momentum 
conservation, it will also have no angular dependence. 

As expected, the dipole harmonic indeed has no sys- 
tematic dependence on trigger particle angle (see Fig. [TJ 
row 1). It turns out that the pt dependence is consistent 
with a contribution from both momentum conservation 
and flow, of roughly the expected size (20|, although a 
more quantitative analysis of this v\ and a proposal for 
how to measure it directly will be presented separately 
25]. 

The most important point here is the lack of angular 
dependence, which is consistent with both effects, what- 
ever their size or dependence on p t . 

It should also be noted that there is no significant con- 
tribution expected from collective flow for the 5th har- 
monic and higher, since the anisotropy coefficients v§ and 
above were found to have a very small hydrodynamic re- 
sponse and very large viscous suppression Q. The ab- 
sence of higher Fourier harmonics provide another piece 
of evidence in favor of the absence of non-flow signals. 



All odd harmonics, then, have no dependence on the 
angle with the event plane. This fact as well as the de- 
pendence on p t is entirely consistent with the properties 
expected from only flow and momentum conservation. 
Any non-flow contribution would have to have these same 
properties to be consistent with data. 

C. Dependence on trigger p t 

The plots in Fig. [1] were all made from the set of data 
where all trigger particles had transverse momentum in 
the range 3-4 GeV, since the statistical uncertainty is 
smaller than the other set of data with trigger p t of 4-6 
GeV, and one can more accurately assess trends in the 
data. Comparing the two data sets, however, provides 
another way to test whether there is a contribution other 
than collective flow. The flow contribution of both sets of 
data should have the same shape as a function of pt, with 
only a difference in normalization from Vn of the trigger 
particle. If there is a non-flow signal that becomes more 
prominent above 4 GeV, it should be apparent. 

Remarkably, the data match expectations from flow. 
Not only do the two sets share the same features as a 
function of angle, but they also have the same pt depen- 
dence. Further, the difference in absolute normalization 
of the harmonics can be used to infer the large pt depen- 
dence of flow coefficients that have not yet been measured 
— comparing normalizations implies that, like t> 2 , v 3 also 
appears to turn over such that the value for particles with 
p t = 4-6 GeV is slightly (~ 10%) lower than at 3-4 GeV. 
This is both an interesting and reasonable result. 

D. Summary of Fourier analysis 

In summary, to be consistent with data, a non-flow sig- 
nal in the analyzed ridge-like dihadron correlation would 
need to have odd Fourier components with no depen- 
dence on the angle <j) s with respect to the event plane, 
a second harmonic with a monotonically decreasing an- 
gular dependence, and a fourth harmonic that decreases 
and then increases with angle. It would also need to 
have a dependence on p t of both the trigger and associ- 
ated particle, in each harmonic, that is the same as that 
expected for flow. In other words, it would have to have 
all the properties of flow. 

A more sensible explanation is that the long-range cor- 
relation data in fact contain measurable contributions 
only from collective flow. 

III. ZYAM BACKGROUND SUBTRACTION 
AND SPURIOUS SIGNALS 

The authors of Ref. [3] came to a different conclusion 
concerning the ridge-like correlation. It is useful to un- 
derstand how they came to this conclusion, and show 
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how the ZYAM subtraction procedure can give mislead- 
ing results. 

As the name implies, the Zero- Yield- At-Minimum pre- 
scription makes specific assumptions about dihadron cor- 
relations, and in a case where such assumptions are er- 
roneous, the subtracted signal provides questionable in- 
formation at best. 

Specifically, one must first assume precise knowledge of 
the flow background — often this is assumed to be just a 
contribution from measured elliptic flow, but sometimes 
this is improved by including quadrangular flow, as in 
Ref. 0- Then, in the presence of non-flow, the overall 
measured correlation is of the form m 



STAR ZYAM 
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nA 



dN 
dA0 



--B pF(A0)] + NF(A0) 



--B 



1 + 2v { 2 a) v i 2 t ' R) cos(2A0) + 2vf ) vf R) cos(4A0) 
+ NF(A^) (7) 



Here, F is the assumed flow background and NF is 
the remaining non-flow signal. The flow background is 
typically assumed to be as in Eq. ([3]), and sometimes in- 
cluding Eq. (j4]). In principle, the input flow coefficients 
(e.g., v n ), should come from an independent measure- 
ment that does not contain a contribution from non-flow 
correlations, but still has the same dependence on flow 
fluctuations. Such a measurement does not exist, and it 
should be noted that it can be a very non-trivial exercise 
to determine the correct input values precisely enough to 
reliably extract a small non-flow correlation (see, e.g., the 
discussion of the vi estimation in Sec. Ill Al) . Also note 
the complete lack of data for odd flow coefficients to date, 
which precludes their use in such subtraction schemes. 

Since the non-flow contribution is unknown, the over- 
all normalization, B, is now a free parameter. To extract 
the desired non-flow signal, one fixes B by assuming the 
non-flow contribution NF is always positive except at one 
or more minima, where it is zero (i.e., zero yield at min- 
imum) : 



NF(A</>) 



dN 
~dA4> 



B[l + 2v { 2 a) v!f' R) cos(2A<£) 

+ 2vi a M t ' H) cos(4A(£)], (8) 



NF(A0 min ) = NF'(A0 min ) = 0, 
\dA(bJ A ^ ia F(Aew)' 



(9) 
(10) 



Since dihadron correlations in, e.g., proton-proton and 
deuteron-gold collisions do typically have close to zero 
yield at their minimum, it may seem reasonable to sus- 
pect that the non-flow contribution to heavy ion collisions 
may also have this form. However, there is no particular 
reason that this should be true. 

More importantly, even if it is true, the ZYAM pre- 
scription can only extract this signal if the background 
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FIG. 2. (Color online) Background subtracted correlation 
■j^—j^ = NF(A0) for trigger p t = 3-4 GeV and asso- 
ciated p t = 1.5-2 GeV. Circles are measured data, while 
solid lines represent a sum of the first four harmonics ex- 
tracted from the unsubtracted data. The left column repre- 
sents ZYAM subtraction with the estimated background from 
Ref. while the right column uses the extracted harmonics 
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V4A- The upper row is for an 



in-plane trigger particle, while the lower row is out-of-plane. 



is known to very good accuracy. Due to event- by- event 
fluctuations, however, this is fraught with difficulty. As 
mentioned, it is not possible to independently measure 
the input flow coefficients v n in such a way that removes 
non-flow correlations but still has the same contribution 
from flow fluctuations. Thus, there will always be model 
dependence and uncertainty in these parameters that, 
small as they may be, compete with the level of preci- 
sion needed to extract the desired signal of interest. If 
any of these coefficients are not correctly estimated with 
enough precision (and certainly if the odd harmonics are 
not considered) the ZYAM procedure will not actually 
"subtract" even the assumed background, and the result 
will have significant contribution from flow. This will be 
true even if the zero-yield-at-minimum assumption itself 
is valid. 

This is illustrated in Fig. [2j The circles in the left col- 
umn represent ZYAM-subtracted data as in Ref. for 
trigger particle in-plane and out-of-plane. The solid lines 
represent the same subtraction procedure, but starting 
with only the first four Fourier harmonics of the unsub- 
tracted data (the right side of Eq. ([1])). Note first that 
all useful information is contained in these first few har- 
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The fact that the double-peak structure on the away 
side depends on <f> s was taken as evidence that triangu- 
lar flow does not represent a major contribution (since 
triangular flow should not depend on event plane orien- 
tation), and by extension that collective flow alone can 
not therefore explain the subtracted signal. As we have 
seen, however, the odd Fourier harmonics indeed have 
no dependence on <f> 3 , just as predicted, and therefore 
they are the same in the upper and lower panels of the 
figure. All dependence on 4> s comes only from the even 
harmonics, which have only been partially subtracted by 
the ZYAM scheme. 

It should be noted that this is not only because vi 
has been underestimated. Consider the right column of 
Fig. [2j Here, instead of using the reported values of V2 
and i>4, I simply use the actual extracted Fourier har- 

(a) (t.R) T r (a) <t,R) t t rpi i 

monies v\ 'v\ ' = V 2 a, v\ v\ = V4A- Ihe sub- 
tracted signal is still not independent of trigger angle, 
even when the actual even harmonics are what have been 
"subtracted". This, of course, is because the effects of 
the even harmonics are still present — the subtraction 
scheme didn't actually remove them, because the ZYAM 
assumption is not true for the remaining odd harmonics. 

A more thorough illustration can be found in Ref. [2|| . 
There, the authors used the AMPT Monte Carlo model 
to investigate the dihadron correlations. By randomizing 
the azimuthal angle of the transverse momenta of initial 
partons in their HIJING initial conditions, or turning 
off jet production completely, they can isolate the effect 
of jet correlations. They indeed conclude that no non- 
flow correlations exist at large relative rapidity in their 
simulations (only effects from "soft hot spots", i.e., the 
initial geometry of the system.) 

More importantly their use of ZYAM subtraction il- 
lustrates several of these problems. They are able to 
independently calculate all v n flow coefficients, and so 
have excellent control over the background used in this 
subtraction (and they include all coefficients, including 
triangular flow). Nevertheless, because of the mentioned 
difficulties, they still find characteristic double-hump sig- 
nals — even when jet correlations are completely turned 
off (see, e.g., the blue dash-dot curve in their Fig. 5 [24j). 

Even in the best of circumstances, then, the ZYAM 
background subtraction procedure can produce results 
that are difficult to analyze. The characteristic shoulder 
structures are seen even with access to the best possi- 
ble knowledge of flow coefficients and essentially no non- 
flow signal (i.e., only very small correlations from reso- 
nance decays), and it is difficult to differentiate between 



this case and when there are interesting non-flow signals 
present. 

A better way to analyze dihadron correlations is actu- 
ally illustrated by STAR in Ref. 0. It could be argued 
that the most reliable way that we have to remove non- 
flow correlations is to introduce a large gap in pseudo- 
rapidity. Therefore, we can simply take the long-range 
correlation (at large A77) as our measurement of flow. 
Then, this flow correlation can be directly subtracted to 
investigate the non-flow signals present at shorter ranges. 
An example of this type of procedure was used in Ref. Q 
to analyze "jet-like" correlations at |A?7| < 0.7. This 
completely removes any assumption about the nature of 
non-flow correlations, except that is negligible at large 
relative pscudorapidity. As we have seen, this is a well- 
supported assumption. 



IV. CONCLUSIONS 

In conclusion, I have shown that long-range two- 
particle correlations are consistent with being entirely 
generated by collective flow. Furthermore, analyzing 
such data with a ZYAM subtraction scheme can be 
misleading. Any ZYAM-based analysis requires precise 
knowledge of the flow background and careful attention 
to effects of flow fluctuations, which show up in every 
Fourier component (not just triangular flow), along with 
a reliance on the zero-yield-at-minimum assumption. I 
argue that it is better to take the long-range correlation 
as our best measurement of flow itself. One can then sub- 
tract this from shorter-range correlations to study non- 
flow effects without any assumption about the nature of 
the non-flow signal — similar to what was done by STAR 
to study jet- like correlations in Ref. Q ■ 
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